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Abstract 

In this paper, we generalize a previous relativistic 1 + 1-dimensional model for two 
mass-less Dirac particles with relativistic contact interactions to the V-particle case. 
Our model is based on the notion of a multi-time wave function which, according 
to Dirac, is the central object in a relativistic multi-particle quantum theory in the 
Schrodinger picture. Consequently, we achieve a manifestly Lorentz invariant formu¬ 
lation on configuration space-time. Our model is constructed to be compatible with 
antisymmetry and probability conservation in a relativistic sense. On the mathe¬ 
matical side, we further develop the method of multi-time characteristics and show 
that uniqueness of solutions follows from probability conservation. We prove that the 
model is interacting and outline how one can understand the interaction as effectively 
given by a (^-potential at equal times. Finally, we answer the question whether Lorentz 
invariant and probability-conserving dynamics can also be obtained when the parti¬ 
cles are confined in a region with a non-zero minimal space-like distance, a question 
relevant for an extension to higher dimensions. 

Keywords: relativistic interactions, multi-time wave functions, zero-range physics, 
Dirac equation, boundary conditions 


1 Introduction 

In a preceding paper [1], it was shown that Lorentz transformations of V-particle wave 
functions in the Schrodinger picture lead to the necessity to consider multi-time wave 
functions, i.e. maps 

: O C X • • • X ^ ‘S, ' — > ifixi, ...jXn), (1) 

N times 

where 5 is a suitable spin space. The crucial point is that the Lorentz transformation of 
a simultaneous space-time configuration ((t,xi),..., (t, xat)) in general yields ((tgxi),..., 
(ty,xy)) where the times are not equal. This reasoning motivates to regard the set 
of space-like configurations as the appropriate domain 11. The multi-time wave function 
generalizes the familiar single-time wave function from Schrodinger’s theory and gives 
back the latter for equal times = t Mk. 
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Following Dirac [2], one usually considers N simultaneous first order partial differential 
equations (PDFs) as evolution equations: 

= Hk^, k = l,...,N. (2) 

Schrodinger’s equation is re-obtained from ([^ by the chain rule for ...,t,XAr), with 

Hamiltonian H = ^^^=1 

The strategy employed in [T] to escape no-go theorems about relativistic interactions such 
as for interaction potentials |3] was to prescribe the free multi-time Dirac equation on 
=5^ and to introduce interaction only by boundary conditions on the set ‘to of coincidence 
points in space-time. This idea is related to the field of zero-range physics |1]. It was 
worked out in detail for two mass-less particles in 1 -|- 1 dimensions (d = 1). In this case, 
an explicit solution of the model was feasible via a generalized version of the method of 
characteristics. This method made it possible to prove existence and uniqueness of classi¬ 
cal solutions. Probability conservation on general space-like hypersurfaces, however, had 
to be checked separately. It led to certain conditions on the conserved tensor current of 
the theory which were equivalent to linear relations between the components of the wave 
function. Together with the requirement of Lorentz invariance, this allowed to formulate 
a general class of boundary conditions. This class was also compatible with antisymmetry 
of the wave function in the case of indistinguishable particles and led to interaction in 
the sense that a generic initial product wave function becomes entangled with the time 
evolution. 

In this paper, we extend the previous model to the A^-particle case. In order to achieve 
a concise formulation, we use a bottom-up approach: instead of starting from all mathe¬ 
matically possible classes of boundary conditions which lead to existence and uniqueness 
of solutions and successively restricting the class of boundary conditions according to the 
various physical requirements, we directly specialize on the case of indistinguishable parti¬ 
cles, extract a class of physically reasonable boundary conditions and then prove existence 
and uniqueness for the resulting model. 

The structure of the paper is as follows. In section]^ we introduce the model as defined 
by its multi-time equations, domain, initial values and boundary conditions. In sec. 
the implications of antisymmetry are discussed and the general solution of the multi-time 
equations is found using a generalized version of the method of characteristics. Sec. |^deals 
with the question how to formulate probability conservation on space-like hyper surfaces. 
A class of boundary conditions is extracted which ensures it. Furthermore, we show that 
the uniqueness of solutions of the multi-time equations follows from probability conserva¬ 
tion. In sec. I^the requirements of Lorentz invariance are checked and the previous class of 
boundary conditions is shown to satisfy them. The main result of the paper is an existence 
and uniqueness theorem (sec. [^. In sec. [^we give a general argument that the model is 
interacting and show that one can regard the interaction as given by an effective (5-potential 
at equal times. Moreover, in sec. [^we answer a different question that was raised in [T]: 
Do consistent Lorentz invariant and probability-conserving dynamics exist on the set Dq of 
space-like configurations with a minimum space-like distance a? This question is relevant 
for the question whether one can also achieve interaction effects by boundary conditions 
in higher space-time dimensions. 
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2 Definition of the model 


Our model is based on a multi-time wave function Q for N mass-less Dirac particles on 
the set of space-like configurations, 

^ := {{ti,zi,...,tN,ZN) ■ {tj-tkf - {zj - Zkf < 0\/j k}, (3) 

in 1 -|- 1-dimensional space-time, with metric g = diag(l, —1). The appropriate spin space 
is 5 = Thus, ijj has 2^ spin components V’i, * = Ij 2^. 

As multi-time evolution equations ([^ we use a system of N mass-less Dirac equations 

...,xn) = 0, k = l,...,N. (4) 

Here, Xk = (tk, Zk), dk,^ = and 7 ^ is the g-th Dirac gamma matrix acting on the spin 
index of the k-th. particle. We choose the following representation: 


7 = 0-1 = 



(T1CJ3 = 



(5) 


where Uj, i = 1, 2, 3 denote the Pauli matrices. This representation diagonalizes Q which 
can be seen by multiplying eq. Q with 7 ^ from the left. This results in 


f d d \ 

\dtk ^^^'^dzkj 




0, A: = l,...,iV, 


( 6 ) 


where a^^k is the third Pauli matrix, CJ 3 = diag(l, —1), acting on the spin index of the k-th 
particle. 

Initial data are prescribed on the set 


T := {(ti, 21 ,..., f^v, zat) G D : ti = • • • = tjv = 0}. 


(7) 


Since has a non-empty boundary one should expect that boundary conditions are 
needed to ensure the uniqueness of solutions. At this point, we leave open the exact nature 
of the boundary conditions. It will be clarified by further considerations about Lorentz 
invariance and probability conservation. 

We summarize the structure of the model as follows: 

{ the system of equations (|^ on D = =5^, 

initial conditions on I, ( 8 ) 

boundary conditions on clD. 


3 Antisymmetry, general solution and multi-time character¬ 
istics 

In this section, we first show how antisymmetry of the wave function for indistinguishable 
particles allows to reduce the domain from N disconnected parts to a single connected one. 
Using a notation for the spin components which is tailor-made for the multi-time equations 
(|^ we show how this facilitates to explicitly determine their general solution. This leads 
to the notion of multi-time characteristics. 
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3.1 Antisymmetry and reduction of the domain 


Following the spirit outlined in the introduction, we make simplihcations wherever phys¬ 
ically reasonable in order to achieve a concise model for which existence and uniqueness 
can be proved elegantly. The hrst simplihcation is the assumption of indistinguishable 
particles. This is natural considering that the particles are not dynamically distinguished 
by eqs. Q alone. Denote the spin components of ijj by ijjsi...sN where each Sj can take the 
values ±1. We write 


V’2 

V’3 


/ iIj - \ 

•0 -+ 

ijj -+_ 


(9) 


V '02^ / \ ^++-++ / 


Indistinguishability implies the following antisymmetry condition for the wave function. 
Let TT G be a permutation. Then 


We now use this condition to relate a solution of (|^ on the different parts of the domain 
Q = ^ (see eq. (§)■ Note that in one spatial dimension, ^ separates into Nl disjoint 
parts which can be classihed according to the relation of the spatial coordinates Zk, e.g. 
Z 2 < zi < < zs < .... Using the permutation group , we write ^ as the disjoint 

union of open sets as follows: 

= u 

where ■= {(U, 21 ,..., ^at, ^at) G-5^ : ^7r(i) < ''' < ^TrjiV)} • (H) 


The crucial point is the following: given a solution of the model, as dehned by ([^ on 
(corresponding to 2:1 < ••• < ^at), antisymmetric continuation via eq. (10) yields a 
solution on provided the boundary and initial conditions are chosen to be compatible 
with antisymmetry. Note that this reduces the possible classes of initial boundary value 
problems (IBVPs) ([^ to an autonomous IBVP on We shall employ this strategy in 
the following. Our new model may be summarized according to (|^ with 5^ replaced by 
^ 1 . 


3.2 Multi-time characteristics and general solution 

Using the notation ([^, we express the diagonalized multi-time Dirac equations (|^ for a 
hxed component 'ipsi...sN follows: 




( 12 ) 


Note that (12) imposes N equations for each of the 2^^ spin components 'ipsi...sN- This 
simple form of the equations allows to hnd the general solution. 


Lemma 3.1 The general solution of eqs. (12) is given by 

...sjv (U) ■2-1 ) • •.) t]\f , Ztv) — /si ...S]v T T 'SW^Al) 

where fsi...si.^- C are C^-functions, si = ±1,..., sn = ±1. 


(13) 
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Proof: The result is obvious if one is familiar with the notation. Simply write out eq. ( 12 ) 
for 'ipsi...sN separately: V’si...sjv = ■■ 






= 0 . 


implies the form (13). 


This 

□ 


The form of the general solution motivates the following definition. 
Definition: Let p = ... t]\f, zn) G . Then we call 


Cfc •— ■2'A: T 


(14) 


the characteristic values at p associated with the component ipsi...sN- 
Furthermore, we define the multi-time characteristics of the components V'si...sjv 

*S'si...sjv(ci, •••jCat) := {{ti,zi,...,tN,ZN) : Zk + Sktk = Ck}. (15) 


With these definitions, one can reformulate lemma 3.1 as follows: the components 
of solutions of (12) are constant on the respective multi-time characteristics (15). Note 
that this implies existence and uniqueness on the domain for an initial value problem 
at ti = ... = = 0, the functions being given by the initial values. However, as 

known from jl], this is in general not true for a domain with boundary such as 


4 Probability conservation 


In this section, we specify an adequate notion of probability conservation on space-like 
hypersurfaces using the conserved tensor current of the multi-time Dirac equations (sec. 
4.1). This enables us to give a sufficient condition on the components of the wave-function 


which leads to probability conservation (sec. 4.2). Furthermore, we prove a general theorem 
showing that probability conservation implies the uniqueness of solutions (sec. 4.3). 


4.1 A geometric formulation of probability conservation 

It is clear that the non-relativistic notion of probability conservation, f d'^xi ■ ■ ■ d'^XN 
|V’P(t, xi,..., t, xtv) = 1 Vt, which heavily draws on a notion of simultaneity, has to be 
generalized in a relativistic context. Building on previous work, such a generalization 
was given and justified in [I] sec. 4]. It makes use of the conserved tensor current of the 
multi-time Dirac equations (Q, i.e. 

jUi-UN V'Tr • ■ with = 0 V/c. (16) 


We define the current form as follows: 

d 

ojj := ^ (-l)^i+-+^^j^i-^^(dx? A...(L<i A---Adxf) 

A • • • A{dx% A ... dx'^ • • • A dxf^) 


(17) 


where (•) denotes omission from the wedge product. The continuity equations (16) im¬ 
ply that the exterior derivative of ojj vanishes, i.e. dojj = 0. The relativistic notion of 
probability conservation on a domain D then reads [J: 


for all pairs of space-like hypersurface^ S, 

^We assume all space-like hypersurfaces to be smooth for the rest of the paper. 


(18) 
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4.2 Boundary conditions from probability conservation 

The formulation via the A^d-form uj together with the property dojj = 0 makes it possible to 
use Stokes’ theorem to extract conditions on ojj and thereby on j which ensure probability 
conservation. 

Lemma 4.1 Probability conservation on in the sense of 

[ ( 19 ) 

for all space-like hypersurfaees S, S' holds if the wave funetion is eompaetly supported 
on all sets of the form S^ H and if 


~ 9 , ( 20 ) 

where 

\= {{tx,zi,... tN, zn) e \3k:tk = 4+i A 2;^ = Zk+i} ■ (21) 


Remark: 1. The assumption of compact support of the wave function (or, alterna¬ 

tively, of suitable drop-off conditions) with respect to spatial directions is needed as 
a technical assumption in the proof. It is reasonable because the multi-time Dirac 
equations have finite propagation speed (see eq. (dU)). Consequently, compactly 
supported initial data imply the desired property. 


2. Note that the wave function is, strictly speaking, not defined on d^i. When using 
values of the wave function at the boundary (such as in eq. (20)), we assume that 
the wave function is continuou^ and refer to the corresponding limit in In this 
way, jumps of the wave function across the boundaries of different are admitted. 
In fact, singularities of this kind are typical for zero-range interactions [H appendix 
J]. 


Proof: We adopt the idea of [H proof of thm. 4.4] and generalize it for N particles. Let 
S, S' be space-like hypersurfaces. We construct a suitable submanifold with boundary in 
order to be able to use Stokes’ theorem. 

Let tjxiz) denote the time coordinate of the unique point p = {tjx{z),z) G S. Let R > 0 
and consider the following set: 


Vr := \{tx,zi,...,tN,ZN)&yi 


3t G [0, l]:\/k:tk = ts(^fc) + t {tT.'{zk) - tT,{zk)) 
and \zk\ < R 

( 22 ) 

Vr is a bounded and closed, thus compact, {N -|- l)-dimensional submanifold of with 
boundary 

OVr = (s^n^i)u((s')^n^i)uMiUM2 (23) 

where M 2 is the subset of Vr with \zk\ = R for some k and 


Ml = VrCS^i. 


(24) 


Because of the first condition in the definition of Vr, a conhguration in Vr is always an 
element of for some space-like hypersurface S. Therefore, it can only be an element of 

^The assumption of continuity is justified in sec. where it is shown that a unique C''' solution exists 
for an appropriate IBVP. 
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Ml C d^i (i.e. light-like) if 3k : = ffc+i and Zk = Zk+i- This implies Mi C 

In the limit i? —)• oo, the integral Uj vanishes becanse of the compact support of the 
wave function. Thus, it follows from the the theorem of Stokes, together with dojj = 0, 
that 


0 = lim / 


doJi = lim 
R-^oo 


OJi = — 


LOn -\- 


OJi + 


OJi 


OVr 






I Ml 


(25) 


The minus sign in front of the hrst integral on the r.h.s. is due to orientation conventions. 
Thus, probability conservation in the sense of eq. (19) holds iff ujj = 0. In order to 
make this integral vanish for ah possible choices of S, S', the condition 


has to be satished. 




(26) 

□ 


Next, we study the implications of condition (26) for the components of the wave function. 
For simplicity, we hrst focus on the special case of equal-time hypersurfaces in a hxed but 
otherwise arbitrary Lorentz frame. 


Lemma 4.2 Let ^i^t ;= {(ti, zi, z^) G = • • • = tiv}. Then the eondition for 

probability conservation on equal-time hypersurfaces St in a particular Lorentz frame, i.e. 
(26) with '^1 replaced by ‘rfi^t, holds if and only if the following condition is satisfied: 

- cr3,k+i) fj{p) = 0 Vp G Mk = 1, ...,N - 1, (27) 


where := {{ti,zi, ...,tN,ZN) G 'ifpt : Zk = zu+i]. 
Furthermore, eq. (27) can be rewritten as 


S 


sk+i\'ipsi...sN\‘^{p) = 0 Vpe^l^J. 


(28) 


(si ,...,sjv)S{— 


Proof: We have to evaluate the condition ^ = 0. Note that for p G 'rfi^t there exists 
a k G {1,...,A' — 1} such that p = {t, zi, ...,t, Zk = z,t, Zk+i = z, ...,t, zn)- 
Next, we calculate according to eq. & recalling that in this case = t and 

Xj = Zj as well as Zk = Zk+i = z. All terms with more than one index p; = 1 in jUi--ui-WN 
vanish because they contain dt Adt = 0. Moreover, the terms with = pfc+i = 0 do not 
contribute, either, as they contain dz A dz = 0. We are left with terms where all indices 
Hj are equal to zero apart from the k-th or the {k -|- l)-th: 


‘-^jip) = —dzi A • • • A dzk-i Adz Adt A dzk+2 A • • • A dzN 

_ ^•o...(Atfe=i)(Aife+i=o)...o^p^ A • • • A dzk-i Adt Adz A dzk+2 A • • • A dz^ 

_ ^^■o...io...o _ ^•o...oi...o^ A • • • A dzk-i Adz Adt A dzk+2 A • • • A dz^ ( 29 ) 


This expression vanishes if and only if the bracket is zero. This yields condition (27): 

0 = (/•■■10...0_^.0...01...0)(^) ^ V’^(p)(a3,fc-CT3,fc+l)V’(p) 

Written out in components, eq. (pOj) reads: 


(30) 


0= X] {-Sklf^si-sJ^ip) + Sk+l\'lpsi...SN\'^{p)) , 

(si,...,siv)e{±i}'^ 


(31) 
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where k was defined above. 

Summands with Sk = Sfc+i cancel out. We are left with 


0 = 


E 


\Hp)- 


^k^^k + 1 


which yields (28). 


(32) 


□ 


We now take the following approach in order to find adequate boundary conditions that 
lead to probability conservation on general space-like hyper surfaces. First, we choose a 
subclass of (28) which turn out to be Lorentz invariant (see the next section) and to 
ensure the existence of a solution (see sec. [^. Then we prescribe the condition on the 
whole set and show that it is indeed sufficient to ensure condition (20) and therefore 
probability conservation on general hypersurfaces. 

It is useful to define the sets 


'^k,k+l ■ — { (^1) 7 ^ 11 3/c . — tfc+l Z}^ — ^k+l } ■ 

One can then write ‘^\ = Uf=i' (see eq. (HI))- 

Lemma 4.3 Let G {~'^i for /c = 1,..., — 1. Then the boundary conditions 


(33) 


/ _L I 

—Sfc+2...SJV ~ ® ¥si...Sfc_i— \-Sk+2---SN 


fk,k+i, k = l,...,N-l (34) 


imply probability conservation on all space-like hypersurfaces in the sense of eq. (19). 


Proof: It was shown in lemma 4.1 that equation (19) follows if 

~ 0. (35) 

We show that this equation indeed holds. Pick a point p € ‘rfi. Then 3k : p £ ^k,k-\-i- 
Condition (34) at this point yields: 


It follows that 


—Sfc+2...Sivl (P) — |V^si...Sfc_i—l-Sfc+2...sjvl (P)- 


jUlU2-.01...UN — jUlU2.-^0...UN (jPj 


(36) 


(37) 


because the expression for the current is diagonal in the components. In the formula for 
Ujip) (eq. 0), we can first sum over the indices pik^Pk+i afterwards over the rest. 
Then there are four possibilities in the summands: 

• {pLk^P-k+i) = (0,0) or (1,1): These do not contribute because the coordinates of the 
A:-th and (fe-l-l)-th particles are equal, say to (t, 2 ), so either dzAdz = 0 or dtAdt = 0 
appears as a factor in the wedge product. 

• (/Ufc, ;Ufc_|_i) = (0,1) or (1,0). One can see that these two factors cancel each other 
because (abbreviating the other factors in the wedge product by A and B) 


jUi-(Pk=o){uk+i=i)-UNji^ /\dtAdzAB + j^^2■■■{uk=l){uk+l=o)...UNJ^ AdzAdtAB 

= {f^■■■^^■■■>^^{p)-f^■■■^°■■■^^^{p))AAdtAdzAB ® 0 . 

Therefore, the probability-conserving property ojj{p) = 0 holds. 


(38) 

□ 
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4.3 Probability conservation implies uniqueness of solutions 


The notion (18) of probability conservation is very powerful. In this section, we work out 
the claim in [1] that is a so-called energy integral and that therefore probability 

conservation implies uniqueness of solutions in a suitable sense. 


Definition: Let S be a space-like hypersurface. We define function spaces 

:= (39) 

Furthermore, we call the solution of the IB VP ^ weakly unique iff for every two solutions 
and every space-like hypersurface S the restrictions of to arguments in 

n P are equal as elements of 'H^\ 


Theorem 4.4 Consider the IBVP (|^ with boundary conditions ensuring probability con¬ 
servation (18) and initial values onX= (So)'^nP, i.e. solution 

is weakly unique. 


Proof: Consider the expression 


ll<^lll := 




U}j{(j)), 


(40) 


where ojj (cp) is the current form constructed from (p according to eqs. (16) and 0- Because 
the Dirac tensor current j is positive-definite and sesquilinear in the wave function, || • ||s 


defines a norm on Ti 


(TV) 


Let 'ip,g) be solutions of the IBVP. Then: ^ S ^ therefore IIV’Is^ ~ 

V^|s„llso = 0- Now let S be an arbitrary space-like hypersurface. Probability conservation 
(18) yields: 

Eo ~ ‘^Ieq 11^0 ~ 0 




Is = 


(41) 


and it follows that as elements of % 


(TV) 


□ 


Remark: The proof of thm. 4.4 suggests that the map 




n 


(TV) 
S' ’ 


V^Ie ^ 


(42) 


which sends the restriction of a solution ip of the IBVP to n P to its restriction to 
(S')^nP, defines a unitary evolution map from one space-like hypersurface to another (see 
also |5l sec. 3]). The analogous view in quantum field theory constitutes the Tomonaga- 
Schwinger picture laizi- 

Note that having in mind a functional-analytic view on time evolution, it might seem 
natural to take the reverse way to define a multi-time evolution, i.e. first defining the 
spaces and a unitary map Vs-s-s'- However, this is not convincing because then 

there may exist S / S' with S n S' / 0 such that 'ip\^[q) / V’|j./(Q') even for g G S n S' 
(see El)- Without additional conditions to enforce 'p^\.^{q) = for g G S n S', 

this would mean that one could not regard the multi-time wave functions and the tensor 
current j as geometrical objects. This, however, may be necessary for a consistent physical 
interpretation (see e.g. m)- 
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5 Lorentz invariance 


The Lorentz invariance of onr model reqnires the invariance of the domain, the mnlti-time 
wave eqnations and the bonndary conditions. Apart from the last point, the invariance is 
already manifest. In this section, we show that also the class (34) of probability-conserving 
bonndary conditions is Lorentz invariant, meaning that the Lorentz-transformed bonndary 
conditions are satisfied as a conseqnence of the old ones. 

The transformation behavionr of spinors nnder a Lorentz transformation K : x ^ x'm the 
proper Lorentz gronp C'_^ is given by 


...,xat) = 5(A) (g) • • • (g) S'(A)';/’(A ^xi,...,A ^xat) 


where 


5(A) = exp (=-^ [7^ 7 I ■ 


(43) 

(44) 


Here, oj is an antisymmetric (1 -|- d) x (1 -|- d) matrix which characterizes A. 

For d = 1, there is only one free parameter /? G M corresponding to a boost in z-direction. 
One obtains 

cosh (3 + sinh f3 0 

0 cosh /3 — sinh /3 


5(A) = 


(45) 


As the matrix is diagonal dne to onr choice of 7 -matrices, it is easy to calcnlate the V-fold 
tensor prodnct in eq. (43). In this way, we find that the components of V' transform as 

N 

(cosh/3 - Sfcsinh/3)V’si..,^j^(A“^xi,..., A'^xat). (46) 

k=l 

This means that one obtains a factor of (cosh/3 — sinh/3) for every pins and a factor of 
(cosh/3 -|-sinh/3) for every minns in the index (si...S 7 v)- Hence, components with an eqnal 
nnmber of pins and minns signs transform in the same way. 


Example: We discnss the case = 3 in order to motivate the general form (34) of the 
bonndary conditions. Consider a bonndary point p = {t, zi,t, z,t, z) G We nse eq. 

(28) to compnte explicitly what the condition of probability conservation amonnts to: 

ujj{p) = 0 77 |V'__+|2(p) - |V'-+-|^(p) + |'0+-+P(p) - IV'+g-PCp) = 0. (47) 


Now we Lorentz transform this condition according to eq. ( |46[ ) nsing the identity 
(cosh /3 — sinh /3) (cosh /3 -)- sinh /3) = 1: 

0 = (cosh/3 - sinh/3) (|V’—+P - \'ip-+-\^) (p')+(cosh,d -F sinh/3) (|V'+-+|^ - |V'++_P) (p), 

(48) 

where y = (A ^{t,zi),A ^{t,z),A ^{t,z)). One can see that this eqnation cannot possibly 
be Lorentz invariant as a whole. Thns, demanding Lorentz invariance, we have to split 
np eq. (47) into two separate conditions relating only components which have the same 
nnmber of pins and minns signs in their indices, i.e. 

!„/. | 2 /„^ |2 


These eqnations are eqnivalent to 


(p) = 0 , 

l^/’+-+P(p) - IV’++-P(p) = 0. 

(49) 

'lp-+-{p), 

V'++-(p) = e*‘^+V'+-+(p)- 

(50) 
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A priori, the phases (p± could be functions of all particle coordinates. However, it is rea¬ 
sonable to demand invariance under Poincare transformations. Then the phases (p± may 
only depend on the Minkowski distances of pairs of particles, (xj — Xj)^{xi — Xj)^. In the 
3-particle case, there is only one such variable, := (ti — — {zi — z'f‘. However, the 

Minkowski distance changes along a multi-time characteristic. However, this would lead 
to a contradiction because the solution has to be constant along the characteristic (see 
lemma 3.1). Thus, Lp± must be constant in order for solutions to exist. A further inves¬ 


tigations of the existence and uniqueness problem shows that even = (/?_ is necessary 
for solutions to exist. 

This yields the following general picture. Exchanging Sfe+i in ijj, 


SI.. ■ SkSk+l-'-^N 


on 


^k,k+i only yields a phase factor which must not depend on the other spin indices (but 
may depend on k). These insights motivate the choice of boundary conditions (34). The 
following lemma shows that this choice is indeed Lorentz invariant. 

Lemma 5.1 The probability conserving boundary conditions of lemma \4.3\ i.e. 

I • (k) 

— Sk+2---SN = e '4’si...sk-i —i-Sfe+ 2 ...siv A; = 1,..., — 1 (51) 


are Lorentz invariant. 


Proof: According to eq. (46), eq. (51) has the same form in every Lorentz frame. Besides, 
the sets ^k,k+i on which the condition is prescribed are Lorentz invariant. □ 

Remark: One may ask if there are other possible choices of boundary conditions which 
lead to ojj = 0 at the boundary and are Lorentz invariant. The example shows that for 
N = 3 we have already found the only one. For A^ > 4 there may exist more complicated 
boundary conditions with the desired properties. However, aiming at a model valid for all 
N > 2, we do not further pursue this question here. 


6 Existence and uniqueness of solutions 

We now come to the main result of the paper: the theorem on the existence and uniqueness 


of solutions for the boundary conditions discussed so far (thm. 6.1). Furthermore, we find 
an explicit formula for the unique solution of a given IB VP within the class defined by eq. 


with H = = 5^1 and boundary conditions (34) 


We start with providing some intuition about the behaviour of the solution for N = 3. 
The main idea is to make use of the fact that the components of the solution have to be 
constant along the multi-time characteristics [ 1 ]. 

Example: For N = 3 the wave function has 2^ = 8 components. According to lemma 


3.1, these are constant along their respective multi-time characteristics. We visualize the 
multi-time characteristics as follows (see fig. Q. One can see from eq. (15) that the multi¬ 
time characteristics are the Cartesian product of V = 3 lines. These lines are plotted 
in the same space-time diagram. Any combination of three points on the different lines 
constitutes an element of the respective multi-time characteristic. The slopes of the various 
lines characterize the associated component iI>si...sn- More precisely, a line with nevative 
(positive) slope for particle k is associated with the appearance of Sk = -|-1 (sfc = — 1 ) in 
the index of 

Fig. a shows a multi-time characteristic S'+++(ci, C 2 , C 3 ) for the component ^/>+++ where 
the Ck are defined by a certain point p = {A, B, C) G (see eq. ®)- '!/>+++ is determined 
at p by initial data because the whole characteristic S+++(ci, C 2 , C 3 ) is contained in .5A\. 
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This can be seen from the fact that every three points on different lines are space-like 
related. Besides, the value V’(p) is determined uniquely by initial data as there exists a 
unique intersection point of 5+++(ci, C 2 , C 3 ) with the surface ti = t 2 = H = 0 , given by 
(0, Cl, 0, C 2 , 0, C 3 ). Thus, we obtain 

'^+++{tl,Zl,t 2 ,Z 2 ,h,Z:i) = 5 r+++( 0 ,Ci, 0 ,C 2 , 0 ,C 3 ), (52) 

where g+++ are the initial values for 



Figure 1: A multi-time characteristic S'+++(ci, C 2 , C 3 ) for the component V’+++- 'S'+++(ci, C 2 , C 3 ) 
is the Cartesian product of three lines which are plotted in the same space-time diagram. Every 
triple of points on different lines, e.g. {A,B,C), is contained in 5^\. 


For a component 'ijjsis 2 S 3 containing plus as well as minus signs in the index, for example 

_the situation is different (see fig. [^. One can see that intersections of the lines 

defining a characteristic _do occur in the diagram. 



Figure 2: A multi-time characteristic _for _depicted for the same configuration as in 

fig-B One cannot trace back the lines to the initial data surface I because one leaves at the 
point P. Instead one has to first make use of the boundary conditions and can trace back the 

multi-time characteristic for ip _which corresponds to the same lines but with particle labels 1 

and 2 exchanged. 


When an intersection occurs, the multi-time characteristic leaves Therefore, a situation 
like in hg. can occur: tracing back the multi-time characteristic to the initial data 

surface, one leaves the domain. Thus, ip-\ _|-(A, B, C) is not defined solely by initial values. 

To obtain the value of ipj^ _we first realize ipj^ _ \-{A,B,C) = V'h_ \-{Pi P^C). 

Then we employ the boundary conditions to obtain ipj^ _|_(A, 5,(7) = ip-\ _ \.{P,P,C) = 

_ y+{P, P, C). The component ip _|-+ is now determined at (P, P, C) by initial data 

in a similar way as before, i.e. ip- ++{P,P,C)=g _^+(0, C2, 0, Cl, 0, C3). Summarizing the 

relations, we obtain: 

1p+-+{tl,Zl,t2-,Z2,t3,Z3) = e*‘^^'^fl'_++( 0 ,C 2 , 0 ,Cl, 0 ,C 3 ). (53) 
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Pictorially speaking, this amounts to exchanging particle labels and picking up a phase 
while leaving the domain on the way back to the initial data surface. 

The considerations above hint at a general idea: it is possible to obtain an explicit formula 
for the solutions of the IB VP (see eq. (56)) by a process of successively tracing back 
components to collisions, using the boundary conditions to switch the component, tracing 
back to the next collision and finally arriving at the initial data. In this way, one can also 
determine values of components with multiple intersections of the lines constituting the 
multi-time characteristic, like in fig. This motivates the theorem below. 



Figure 3: A multi-time characteristic with several intersections for the component ■0^_ 


Theorem 6.1 Let m G N and choose initial data gj G C™ (Xn.5^i,C) Vj = 
such that they also satisfy the boundary conditions, i.e. 

—Sfe+2...sjv = —hS)c+2---«iV 011 21 n k = l,...,N—l (54) 


and let this transition be C'™'. 

Then there exists a unique solution if G C'™' (C^)®'^) of the IBVP 

boundary conditions (34) and initial data 


with LI = 5T\, 


lfsi...SNiO,Zl,...,0,ZN) = gsi...SNizi,---,SN), Zl<---<ZN- 


(55) 


If all characteristic values Ck = Zk + s^tk are different, the components of the solution are 
explicitly given by 

Z'X, ^“-itN ■} ZN^ e l ^ •••; ; (^^) 

where vr is the permutation with 


VS1...SJV 

below. 


0-77(1) V * * * (51^) 

is the phase which is uniquely determined by the of eq. (34) via the definition 


If some of the Ck are equal, if, 


S1...SN 


is given by the continuation of eq. (56). 


In addition, the model ensures probability conservation on general space-like hypersurfaces 
and is Lorentz invariant. 


Definition: A pair [k, 1) G {1, ..., with k < I is said to be a collision of a transposition 
TT G Sn iff 7r{k) > tt{1). 


Definition: The phases 


^si...sjv 


appearing in thm. 6.1 are given by the conditions: 


1 . 


^id 


= 0 , 
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2. Let Tfc be the transposition of k and k + 1. If vr can be decomposed as vr = o fj 
where cr is a permutation with fewer collisions than vr, then 




Sl...Sk + lSk-"SN 


+ 


(k) 


(58) 


Lemma 6.2 The phases ({)[ 




exist and are uniquely determined. 


Proof: We proceed via induction over the number of collisions. 

Induction start: If tt has no collision, vr = id. Thus, the phase is determined 

uniquely by the first condition in the definition. If vr has exactly one collision, then it is 
just a transposition of neighbouring elements, so there exists some k with vr = r*, = o id 
and the phase is uniquely determined by (58) as ~ ■ 

Induction step: Assume that all phases permutations with n > 1 collisions 

are uniquely determined and let vr have n + 1 collisions. It is known from the general 
theory of permutations that there exists at least one permutation a with n collisions and 
a neighbouring transposition Ts such that vr = Ts o a. However, it may be possible to 
decompose vr in two different ways: 


vr = Ts o a = TfcOK (59) 

where s,k € {1 ,... ,N}, s ^ k and a, k are permutations with at least n collisions. In 
order for these two permutations to have one collision less than vr, we see that (fc, A: + 1) 
and (s, s + 1) must be collisions of vr. 

To show that despite the different ways of decomposition, the corresponding phases are 
uniquely defined, we make use of the fact that the phases need only be defined for a cer¬ 
tain type of permutation. To characterize them, we prove an auxiliary claim: in the above 
situation, Tg commutes with because |s — A:| ^ 1. 


Claim: Let ... ,t]sr, zjsf) G such that there is a collision, i.e. a pair 

(o, b) with a < b and Ca > Cb. Then one of the following two possibilities holds: 


J either Sq = -|-1 A = — 1 A ta > 0 A 4 > 0 
\ or Sq = —1 A Sfo = -1-1 A ta < 0 A 4 < 0 


(60) 


Proof of the Claim: We know that a < b, Ca > Cb and Za < Zb- We show that Sa = +1 
implies Sb = — 1 A ta > 0 A > 0; the second case follows analogously. 

If Sa = +1, then 


Ca ^ Cft 

4A Za T ta A* Zb ~\~ Sbtb 

ta (^btb) ^ Zb Za — \zb Za\‘ (01) 

If now Sb = +1, this would be a contradiction to the points {ta,Za) and {tb,Zb) being 
space-like separated. Hence, = — 1, so we have: 

\zb-Za\ < ta + tb = \ta + tb\. (62) 

This implies that ta and tb cannot both be negative. So assume that one of them is 
negative, w.l.o.g. ta > 0,4 < 0. But then |ta — tb| > |ta -|- tb| and 

\Zk-Za\ < |ta+4| < |ta-tb|, (63) 
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which also is a contradiction to the points being space-like. Thus, one must have 
ta > 0, tfe > 0, which proves the claim. 


Because of the specific sign combinations that allow for collisions, the claim shows that if 
(s, s -|- 1) is a collision, neither (s — 1, s) nor (s -|- 1, s -|- 2) can be one. Therefore, | A: — s| > 2 
and TfcjTs commute. 

We use the commutability of and Tg to define a third permutation 

p := TkOTgOTT = Tfc o fj = TgO K, (64) 


which by construction has n — 1 collisions, i.e. one less than a and k. This means the 
seemingly different representations of , 




-|- and 

^Sl...s^ ~ V"Sl...Sfc_HSfc...SAr ' (65) 

are in fact equal. This can be seen from the fact that the different ways of decomposing vr 
via eq. (59) yield (using 


Sl...Ss + lSs...Sjv 


4’si...Sff 4’si...Ss + lSs...Sff + ^ ^Sl...S,, + iS^....Ss + lSs...Sff + ^ 

= ( 66 ) 

This finishes the proof of uniqueness of the phases because by the induction hypothesis, 
the phases associated with p, a and k exist and are uniquely determined. □ 


Proof of the theorem: The points of Lorentz invariance and probability conservation 


are clear from lemma 5.1 and lemma [4~^ respectively. Furthermore, we already know that 
the uniqueness of solutions in a weak sense follows from probability conservation by virtue 
of thm. 4.4 If the function defined by eq. (56) is indeed m times continuously differentiable, 
it follows from continuity that it is also unique as a C'™'-function. 

Thus, it only remains to show that the function given by (56) is indeed a classical solution 
of the IBVP. In order to prove this, the following four points have to be verified: 


1. Differentiability: We need to prove that ip G (j^i, (C^)®^). As the initial values 
satisfy gj G C™ (X, C) Vj = 1,..., 2^, this property is inherited by ipj via the charac¬ 
teristics. To see this, note that eq. ( |56[ ) just makes use of a translation of the initial 
values along straight lines in the {tk,Zk) spaces. 

However, we need to consider those points separately where the permutation vr 
changes. This exactly happens when at least two of the characteristic values cj 
are equal. But then the C^-property of ip is assured by the requirement that the 
initial values must satisfy the boundary conditions (eq. (|54[)) and that the transition 
shall be C™. 


2. The function defined by eq. (56) solves the system of Dirac equations in This 


follows from lemma |3.1| because the components of the solution are indeed constant 
along the respective multi-time characteristics and only depend on the characteristic 
values Cfc. 


3. The initial conditions (55) are satisfied: At a point (0, zi, 0, 2 : 2 ,..., 0, zat) GXn=5^i 


we have Ck = Zk VA: and thus ^ ••• ^ c^(w) is fulfilled for vr = id. 

Therefore, formula (56) reduces to 




S1...SJV 


( 0 , 2 ;i,..., 0 , 2 ;Ar) = (ci,..., cat) 


(67) 


which is equivalent to (55). 
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4. The boundary conditions (34) are satisfied: Let k G — 1} and (ti, zi, ...,tk = 

t,Zk = z,tk+i = = z, ...,tN, zn) G 'Tfk,k+i- We consider two components of 

where only the /c-th and [k + l)-th sign is exchanged, or more formally: Let 
(si,... sat), (si,... sat) G {-1, +1}^ with Si = s; V/ ^ {k, k + l} and (sfc, Sk+i) = 
(+, —) = (sfc+i, Sfc). The respective characteristic values are given by Ck = Zk + Sktk 
and c-k — Zk T Sktk- 

Now observe the property ci = ci'il ^ {k,k + 1} and Ck = Cfc+i, Cfc+i = Ck- Let vr 
be the permutation that leads to Cjr(i) < • • • < c^(Af)- The permutation a needed to 
achieve 5^(1) < • • • < Co.(Ar) is given by cr = o tt, and it has one collision less than 
vr with respect to the indices Sk- Inserting this into eq. (56) yields 


YSI...SN — ^ 


^SI...SN — c. - 






l|56ll _i^(k) 


— e 


ih 


( 68 ) 


This shows that (34) is valid. 


These four points establish existence; the function given in (56) is indeed the solution of 
the IBVP. □ 


Remark: 1. Uniqueness of solutions can also be proven differently than by invoking 

thm. |4.4[ namely by directly showing that every solution of the IBVP has to fulfil 
eq. (56). A possible proof is via induction over the number of collisions. 


2. Note that on purely dimensional grounds it is remarkable that solutions for the 


IBVP with boundary conditions (34) do exist. Because the dimension of ‘^k,k+i is 
{N — I)(I + d) which for N > 1 + d is greater than Nd, the dimension of the initial 
data surface X n one might have suspected this not to be the case. 


7 Interaction and effective potential 


In addition to the mathematical and physical features already established, we now prove 
in this section that our model is interacting. Moreover, we outline how the interaction can 
be described effectively at equal times using self-adjoint extensions of the free two-particle 
Dirac Hamiltonian. 

In |1], a general criterion for interaction was given. We call a physical model interacting 
iff it generates entanglement, i.e. if there exist wave functions that are initially product 
states and become entangled during time evolution. Note that for the antisymmetrized 
wave functions we are considering, a product means “wedge product”. We now present a 
simple argument why our model is interacting in this sense. 


Lemma 7.1 The model defined by ([^ with H and boundary conditions (34) is 

interacting if there exists k G {1 ,... ,N — 1} with vr. 


Proof: W.l.o.g. k = 1. Let the initial conditions be such that is a product wave 
function. In particular, this means that there exist functions a,l3,'y,6 G C) and 

C G C™(M^-2,C) with 

9+-+-+{zi,-,zn) = a{zi)l3{z2)C{z3,-,ZN) 

9-++-+{zi,-,zn) = j{zi)6{z2)C{z3,...,zn) (69) 
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for 2 i < • • • < zn- 


Antisymmetry (10) implies 


a{zi)l3{z2) = -'y{z2)S{zi). 


(70) 


Consider the solntion at a point p = {t,zi, ...,t,Z]\f) G with common time t > 0. The 
anxiliary claim in the proof of lemma | 6 . 2 | implies that the characteristic valnes at p with 

respect to the component 'ijj^ _|_are in ascending order iS zi + t < Z 2 — t. Thns we 

can nse formnla (56) to obtain _|_ ^{p), with the permntation tt being the identity if 


zi < Z 2 — and the transposition n if zi > Z 2 — 2t. Written via the Heaviside function 
0 , this yields 

'0+_+...+ (p) = c/+_+...+ (ci,...CAr)0(z2-zi -2t) 

«(i) 


+ 5-++-+(c2,Ci,C3,...,CAr)0(2t+ Zi - Z 2 ) 

= a{ci)l3{c2)C{c3, cat) 0 {Z2 - Zi- 2t) 

+ e*‘^''^7(c2)(5(ci)C(c3, cn) 0 (2f + zi - Z 2 ) • 


(71) 


Using (70), the expression simplifies to 


V'+-+-+(p) = a(ci)/3(c2)C(c3, • • • ,CAr) (^0(c2 - Cl) - e*‘^'^^0(ci - C 2 )) . (72) 

This function contains the Heaviside function of a combination of t, zi and Z 2 in a non- 
factorizable way. The 0-function cannot be left away for general initial values (as might 
be the case if they were zero in some regions). Furthermore, because of the prefactor 
of the second summand, we cannot write it as a product as long as D 


Effective single-time model: In the following, we show how an effective single-time 
model can be obtained from our model when considered at equal times = ^2 = Even 
though manifest Lorentz invariance is lost for a single-time model, we consider it instructive 
to connect with this familiar setting. 

For simplicity, let N = 2. We denote the single-time wave function by xizi,Z 2 ,t) := 
zi,t, Z 2 )- Then the single-time model is given by the domain {{zi,Z 2 ,t) G : 2:1 7 ^ 


Z 2 }, initial data at t = 0, boundary conditions (34) (with ti,t 2 replaced by t in all the 
constructions) and the wave equation 


■ dx 


d d 

cr3<8)lw-hl<8)0-3— )x = Hx- 

(JZ\ OZ 2 


(73) 


Note that eq. (73) is obtained from the multi-time equations (|^ by the chain rule. 

We introduce new coordinates u = ^(zi—Z 2 ) and v = ^{zi+Z 2 )- The Hamiltonian becomes 


H = -idiag{d^,du,-du,-dy). 


(74) 


The boundary condition (34) can be reformulated using antisymmetry, i.e. X 2 {u,v,t) = 


—X 3 {—u,v,t), with the result 


\im X 2 iu,v,t) = lim-e ^‘^X 2 {u,v,t), tG 

n/'O •u\0 

limx3(u,n,f) = lim-e"^X3iu,v,t), t£] 

n/'O u\0 


(75) 


The components xi and X 4 evolve freely and have to be continuous and zero at n = 0 


because of antisymmetry. For p = tt, eq. (75) also reduces to the condition of continuity. 
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In that case, the model becomes free - in agreement with lemma 7.1 


The bonndary conditions (75) can be implemented in a rigorons fnnctional-analytic way 

Moreover, the nnitary 


by a family of self-adjoint extensions of H parametrized by (/? |in) 
gronps generated by these self-adjoint extensions can also be obtained as the limits of 
nnitary gronps generated by Hamiltonians Hn = H + Vn with 


Vn{u) = diag(0,-K(«),0), (76) 

where the potentials converge to the (5-function, Vn{u) —)■ (tt — ip)6(u), in the distributional 
sense. In essence, this is what the physics literature shows m I12|R Recalling u = ^{zi — 
Z 2 ), our model can therefore be considered a relativistic version of the 1 -|- 1-dimensional 
multi-particle Dirac equation with a spin-dependent S(zi — 2;2)-potential for ip ^ n. 


8 Non-existence of solutions for configurations with a mini¬ 
mal space-like distance 


In this section, we leave the setting of the model (|^, just retaining the multi-time Dirac 
equation (Q, and focus on a different aspect: the question whether a consistent Lorentz- 
invariant and probability-conserving dynamics exists on the domain of space-like con¬ 
figurations with a minimum space-like distance a. 

The relevance of the question is motivated as follows: one may ask whether an immediate 
generalization of the model (|^ to higher dimensions is possible which still produces inter¬ 
action. This question was answered negatively in |1], the reason being that for d > 1 the 
probability flux across the boundary vanishes without the need for boundary conditions 
and thus, by theorem |4.4[ uniqueness of solutions immediately holds. The model then 
corresponds to the free case. The main point is that for d > 1, the set of coincidence 
points in space-time is too low-dimensional to have impact on the dynamics. Therefore, it 
suggests itself to ask whether a change in the domain (and thereby of its boundary) can 
be made such that the dimension of the set across which the probability flux could leave 
the boundary is increased. A natural choice is the set of a-space-like configurations (here 
for N = 2): 

= |(ti,xi,t2,X2) G X : (ti - ^2)^ - (xi - X2)^ <-0:^1 . (77) 

dS^a has dimension 2d -|- 1 and its intersection with S x S, the set appearing in the proof 
of probability conservation |T], has dimension 2d which is sufficient to have impact on the 
dynamics (the reason being that Wj is a 2d-form). Compare with the dimension 1 -|- d of 
the set of coincidence points. However, because itself has a dimension greater than 
2d, the dimension of the initial data surface, the question arises if there is a consistent 
dynamics on it at all. In the following, we approach this question for the simplest case, 
d = 1, for which the necessary mathematical tools are available, and show that the answer 
is negative. 

First we show that there can only be one kind of boundary conditions with the desired 
properties. In a second step we then prove that the corresponding IB VP on 5^^ does not 
possess non-trivial solutions. We make use of the following definitions: 


= {(t,Zi,t,Z2) G y’a '■ - ^2 > 0} , 

= {(t,Zi,t,Z2) G '■ - Z2 < 0} . (78) 


®Note that |11[ I12| 
theless possible since, 


consider a one-particle Dirac equation. A comparison with these papers is never- 
as evident from eqs. (741, 


(761 the single-time equation (731 with additional po¬ 
tential (761 decouples and x := (XZiXs) satisfies a one-particle Dirac equation in u, i.e. idtx{u,v) = 
{-iu^du G aiVn{u)]x{u,v). 
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We have: U 


Lemma 8.1 Let a > 0 and N = 2. For the multi-time Dirae equations Q on the domain 
there exist no other Poineare invariant boundary eonditions whieh lead to probability 
conservation on every spaee-like hypersurface and which are compatible with antisymmetry, 
than the ones given by: 

V’+-(P) = e=^*‘^'0-+(p) (79) 


with a fixed ip G (—vr, tt]. 

Proof: Let p = {t\, Zi,t 2 -, Z 2 ) G dS^a- Because the two points {ti,zi) and {t 2 -,Z 2 ) are 
space-like separated, there is a Lorentz frame with ti = t 2 - We work in this frame, so we 
can write either p = {tp,z,tp,z a) or p = (tp,z,tp,z — a). The idea is to use Stokes’ 
theorem in a similar way as in the proof of lemma [TT] to obtain a condition for probability 
conservation on equal-time hypersurfaces 'Eri,'^T 2 the considered Lorentz frame. Here, 
w.l.o.g. Ti < T 2 . Let 


V := {{t,zi,t,Z 2 ) e y'a ■ n <t < T 2 } 


(80) 


V plays the same role as Vr in eq. 


for R sufficiently large. Following the strategy of 


the proof of lemma 4.1, one deduces 


0 = dujj= ujj 
JV JdV 


(81) 


Note that in contrast to the proof of lemma 4.1 but similar to jl] proof of thm. 4.4], 
there now exist two connected components of the domain Therefore, probability 

conservation in the form 


LOi = 


UJ: 






is equivalent to 


UJn = 


OJi 


(82) 


(83) 


JMW JmC^) 

where = {{t, zi,t, Z 2 ) G 8^01 '■ maxjzi, Z 2 } = zj A ti < t < T 2 } for j = 1, 2. Observe 
that from (t, zi,t, Z 2 ) G it follows that zj = z^-j a. Furthermore, antisymmetry 
implies: 

ujj{t, z,t, z + a) = —ujj{t,z-\-a,t,z). (84) 

This can be seen from the fact that on ^1, Uj = (Itf _tP ~ IV’H_P) dt A dz (see the proof 

of lemma 4.2). 


Inserting (84) into eq. (83) allows us to conclude: 

L 


ujj = 0. 


(85) 


Jmw JmA) 

As this relation must hold for every Ti,r2, we must have ojj{p) = 0. In components: 

IV’+-(p)l^ - IV’-+(p)P = 0 fi+-{p) = (86) 

where (p : dS^a (—vr, vr] could in principle be a function which is not constant. 

Because p is an arbitrary boundary point, this equation must hold on the whole of dd/’a- 
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Moreover, the requirement of Poincare invariance has the consequence that (p{p) has to be 
locally constant (see the example preceding lemma (4.3)). The domain 5^^ has the two 
connected components 5^^ and by antisymmetry one obtains: 




Thus, indeed no other boundary conditions than (79) are permitted. 


(87) 

□ 


Remark: 


1. A similar proof for distinguishable particles shows that in this case another 


possibility appears: the two contributions in eq. (83) could cancel instead of vanishing 
individually. However, this cancelling is not physically sensible because it would 
imply a non-vanishing current from to 5^^ and vice versa. Provided the Born 
rule holds, the particles could then swap place instantaneously. 


2. The boundary conditions (79) do indeed imply Poincare invariance and probability 
conservation. However, this will not be shown explicitly as they do not lead to the 
existence of dynamics (see the following lemma). 


Lemma 8.2 Let a > 0 and consider the IBVP given by 

i'ykdk,iJ.il^{ti,zi,t 2 ,Z 2 ) = 0 for A: = l,2, 

^/^(0,zi,0,Z2) = g{zi,Z2), 

■0+- = _I- on dy’a 


\zi - 2 : 2 ! > a} 


( 88 ) 

is 


on the domain 5^^. Here, (p G (—7r,7r] and g : {(zi,Z 2 ) G 
supposed to be a -function. 

Then, if there exist real numbers ai < 5i < 02 < 62 with g^ (ai,a2) 0 _(^1,^2) or 

g |_(ai,a2) 0 g ^(61,62) ihe IBVP (88) does not have any -solution. 

Proof: Assume that there exist real numbers ai < bi < 02 < b2 with g _ (01,02) 0 

g _ |_(6 i, 62)- The case of g^ _is similar and will not be shown explicitly. Suppose that if 

is a solution of (88). We obtain a contradiction by constructing points (ti,?/i,^2;2/2) and 
(si, xi, S2) 3^2) £ which lie on the same multi-time characteristic with respect to the 

component if^ _(see fig. Q. 

The construction proceeds as follows: 

1. Choose a point {ti,yi,t 2 , 2/2) on the same multi-time characteristic of if^ _as (0, oi, 0, 02 ) 

and on the boundary of .50^, i.e. 



= 221-^1 
= 222+^2 
= (221 - 222)^ - 

This in particular implies: 

0 -+(H, 22102,222) = g-+{ai,a2). 


(89) 


(90) 


2. Consider the set of points (si, xi, S2, X2) on the same multi-time characteristic as 
(0, 61,0, 62) and on the boundary of i.e. 


This means 


61 = Xl — Si 

62 = X2-[- S2 

(si - 32 ^ = (xi - X2)^ - a^. 

0 _+(si,xi, S2,X2) = g-+{hi,b2). 


(91) 


(92) 
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3. Now select a point (si, xi, S 2 , X 2 ) on the same multi-time characteristic as (ti, yi, ^2; 2/2) 
with respect to the component _, i.e. 


f xi-hsi = yi + h, 
\ X2 - S2 = 2/2 - t 2 . 


(93) 


This implies that the value at (ti,2/I5^2;2/2) can be obtained in two different ways: 
firstly by using the boundary condition at that point and secondly by going along the 
characteristic surfac^to (si, xi, S2, X2) and using the value from there. In formulas: 


ip+-{ti,yi,t 2 ,y 2 ) 221,^2,2/2) 
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(94) 


V'+-(ilj 221,^2, 2 / 2 ) =■ 'lp+-{si,Xi,S2,X2) 

e*‘^V'_+(si,Xi,S2,X2) 

® 9 -+{bi,b2)- (95) 


Thus: 


9-+{bi,b2) 

in contradiction to the assumption. 


5'_+(ai,a2), 


(96) 



Figure 4: Construction in the proof for values ai = 1, 61 = 2, 02 = 5, 62 = 6 and a = VQ. 
The points are Aj = (aj,0), Bj — (bj,0), and Yj = {yj,tj), Xj = {xj,Sj) for j = 1,2. The 
black hyperbola consists of points with space-like distance a to Yi and the grey one of thoses 
with space-like distance of a to Xi. The configurations {Xi,X 2 ) and (15,12) he on the same 
multi-time-characteristic, comprised of the Cartesian product of the two solid black lines. 


This proves the claim, provided the points we use do exist. Indeed, the combination of the 
eight equations (89), (91) and (93) with eight unknowns leads to rather lengthy quadratic 
equations the general solution of which can be found in appendix]^ One explicit solution 
is given in the figure. □ 


Remark: The lemma shows that the most general Lorentz invariant and probability- 
conserving IB VP (88) on is over-determined. Eq. (96) shows that the only admissible 


^One may wonder how it is possible to have a path connecting the two points which neither leaves 
the characteristic nor the domain. This is achieved as follows. Concatenate the two linear paths from 
(ii, yiT2,2/2) to (ti, yi,S 2 ,X 2 ) and from (2i, 1/1, S2,0:2) to (si, *1, S2, X2), so first move the right point from 
I2 to X 2 and afterwards the left from 1! to Xi. One can see from the hyperbolas in figure]^ that this 
path only leaves 5^a at its endpoints. 
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initial data are those for which g _^ is constant (and thus also _). Due to normalization, 

this constant has to be zero. The two other components are exactly those which are 
not affected by boundary conditions. Moreover, it becomes clear from the proof that the 
problem originates from the too high dimension of which implies (regardless of initial 
conditions) that certain components of the wave function have to be constant on sets like 
the initial data surface. One cannot avoid this problem by simply prescribing boundary 
conditions only on a part of the boundary due to the requirement of Lorentz invariance. 


9 Discussion 


In this work, we have developed a rigorous, interacting and explicitly solvable relativistic 
multi-time model for N Dirac particles in 1 -|- 1 dimensions. The main results are (a) the 


extraction of the class (34) of boundary conditions which are compatible with the require¬ 


ments of antisymmetry, manifest Lorentz invariance and probability conservation, (b) the 
proof that uniqueness of solutions of the multi-time equations follows from probability con¬ 
servation on space-like hyper surfaces, as well as (c) the proof of the existence of dynamics 
and the explicit formula for solutions. Concerning (a), we believe that this class is the only 
one compatible with the physical requirements which can be formulated for general N (see 
the remark at the end of sec. [^. Furthermore, we showed that the interaction by boundary 
conditions on sets where the space-time coordinates of two particles coincide can, at equal 
times, be effectively regarded as given by a spin-dependent 5-potential. 

Our results show that even the strictest requirements of relativistic invariance, as em¬ 
bodied by the multi-time formalism, can be rigorously satisfied. This makes clear one 
more time that direct relativistic quantum-mechanical interactions are not generally im¬ 
possible. No-go theorems about relativistic interactions such as (Sj [13] rather rule out 
only specific mechanisms for interactions. The possibility to explicitly solve the model is 
instructive for illustrating how the various physical requirements can be met in the multi¬ 
time formalism. The model therefore also has a certain pedagogical value. Moreover, it 
might serve to explicitly test general claims about relativistic quantum mechanics. The 
fact that the model possesses a conserved tensor current furthermore ensures compatibility 
with realistic quantum theories such as relativistic GRW models nmsi and relativistic 
Bohmian mechanics ElEl I16j . This is of particular interest because these theories have 
so far only been formulated for the non-interacting case, which is, as our model shows, 
not due to a problem inherent in these theories but only due to the fact that no rigorous 
interacting relativistic multi-time theory existed before. 

Motivated by the question whether a generalization of the model to higher dimensions can 
be achieved via the introduction of a minimal space-like distance between the particles, we 
first treated the question of the existence of dynamics for these a-spacelike configurations 
for d = 1. The result was negative. This leads us to believe that also in higher dimen¬ 
sions a consistent dynamics on the domain of a-spacelike configurations does not exist. 
See, however, HZI for an analysis of a different multi-time model for two interacting Dirac 
particles in 1 -|- 3 dimensions. 

Concerning further generalizations of our model, it should be possible to include non-zero 
masses. However, one then looses the possibility to explicitly solve the model and conse¬ 
quently a change of the mathematical tools is required (see the discussion of [I] ). A further 
hint that the inclusion of masses should be unproblematic is given by the connection with 
5-interactions outlined in sec. Once a self-adjoint extension has been found which im¬ 
plements the 5-interaction, one can always add a bounded symmetric operator such as a 
mass term to the Hamiltonian. A different generalization would be the case of variable 
particle numbers. Similarly to Fock space, one could then consider the model on the dif- 
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ferent A^-particle sectors and trying to relate them (see for a discussion of multi-time 
quantum field theories). If successful, this would yield a rigorously interacting multi-time 
QFT model in 1 -|- 1 dimensions which, in addition, might be explicitly solvable. 
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A Explicit formulas for the points used in the proof of lemma 
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In the following we give the solutions of the eight equations (89), (91) and (93) which are 


1 


—ai -|- 6i -|- 2(^2 ~ 25i -|- 62) ~ 2^’ 


1 


1 


1 


1 


used in the proof of lemma 8.2 

yi = «i + ^ {-ai + hi + ^(02 - 26i + ?i2) - , 

2 ti 

y2 = “1 + ;^ ( -ai + + ^(02 - 26 i + b2) + ^^ ] , 

t 2 

Xl 

Si 


X 2 = 


S2 = 


02-62 + 2 (a^ -bl + 26162 - 62 + (62 - 61) (2(02 - 26 i + 62) - 2?)) 
" ^ ( 46 i - 462 + 2(02 - 26 i + 62) - 20 ^ ’ 

61 + ^(02 - 26 i + 62) + 

-(02 — 26 i + 62) + 

62 - - 6^ + 26162 - 6| + (62 - 61) (2(02 - 26 i + 62) - 2'?) 

(261 — 262 + (02 — 26 i + 62) + ’ 

- 6f + 26162 - 6| + (62 - 61) (2(02 - 26 i + 62) - 2^) 

(261 — 262 + (02 — 261 + 62) + 


where 


C = 


(62 - 02)^(6i - oi) + 4a2(52 - 02) 


oi - Ol 

The radicand is positive since oi < 61 and 02 + 62. 


(97) 


(98) 
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